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ABSTRACT 

We present a multi-wavlength infrared study of the nearby, edge-on, spiral galaxy 
' NGC 891. We have examined 20 independent, spatially resolved IR images of this 

\ galaxy, 14 of which are newly reduced and/or previously unpublished images. These 

images span a wavelength regime from A 1.2 /im in which the emission is dominated by 
cool stars, through the MIR, in which emission is dominated by PAHs, to A 850 /im, 
in which emission is dominated by cold dust in thermal equilibrium with the radiation 
field. The changing morphology of the galaxy with wavelength illustrates the changing 
dominant components. We detect extra-planar dust emission in this galaxy, consistent 
with previously published results, but now show that PAH emission is also in the 
halo, to a vertical distance of z ^ 2.5 kpc. We compare the vertical extents of various 
components and find that the PAHs (from A 7.7 and 8 /im data) and warm dust (A 24 
/im) extend to smaller z heights than the cool dust (A 450 /im). For six locations in the 
galaxy for which the S/N was sufficient, we present SEDs of the IR emission, including 
two in the halo — the first time a halo SED in an external galaxy has been presented. 
We have modeled these SEDs and find that the PAH fraction, /pah, is similar to 
Galactic values (within a factor of two), with the lowest value at the galaxy's center, 
consistent with independent results of other galaxies. In the halo environment, the 
fraction of dust exposed to a colder radiation field, f co id, is of order unity, consistent 
<^> with an environment in which there is no star formation. The source of excitation is 

ON . likely from photons escaping from the disk. 
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1 INTRODUCTION large molecules can be present at such high altitudes above 

. . , the plane and also their source of excitation. Presumably, 

In a wide variety of galactic environments, polycyclic aro- „ TT , . ,_, TT , . . . .. , . 

. , ° ,„ . . . lar-llV (lUV) emission from m-disk star forming regions 

matic hydrocarbons (PAHs) are the dominant emitters m . , ., ,. TT c . ot-^ • n 

/ i provides the excitation. However, a fit to a SED is really 

the Mid-Infrared (MIR) from about A 5 - 15 um (Vogler et , . , , , ,, , . . , . r 

, ' _ „ , , , V f required to understand both the contributions from vari- 

al. 2005, Smith et al. 2007, Galliano et al. 2008b). At longer . . ^ , . , TT , . , 

ous emission components (cool stars, PAHs, classical grains, 

wavelengths approaching the Far-Infrared (FIR), thermal . N „ . „ ,, , . . , „ 

° ° ., , nm . etc.) as well as the stellar component that is required for 

emission dominates the spectral energy distribution (SED) ,, .. „ , on-r^ n , , , i 

. . the excitation. Better yet, a SED fit at different focations m 

(Galliano et al. 2008a . PAHs have been observed m many , . „ . , , . .. , . , 

, ' a galaxy, especially at different locations m the disk and/or 

galaxies m locations often associated with star forming re- , , r , , , , ' , 

. ° m extrapfanar gas, can further hefp us to understand the 

gions. However, new discoveries of PAH emission extending , . , , , , . TT -. ■ . , n 

compfex interpfay between various IR-emittmg mtersteffar 



medium (ISM) components and the origin of the excitation. 



to significant distances (6 - 10 kpc) from the disks of several 
galaxies (Irwin & Madden 2006; Engelbracht et al. 2006; Ir- 
win et al. 2007) raise important questions as to how these We have therefore chosen to study MIR-submm emis- 
sion in a very well known edge-on galaxy, NGC 891 (Fig- 
ure 1), that has known extraplanar emission in a variety 

i T ,,. , ... „„,„,, , j «t I, . of wavebands (e.g. H I, Ha, dust and CO - See Sect. 1.1). 

in this paper, we adopt the term PAH to describe the carriers v ° ' ' ' 

of MIR spectral feature emission, whereas we refer to thermal Since there is so much spatially resolved IR data available 

continuum emission as dust emission. Note that some authors f or this galaxy, it is an excellent target for modeling its SED 

consider PAHs to be low mass carbon grains in their treatment at different positions in both the disk and halo regions. For 

of interstellar dust (e.g. Draine & Li 2007b). this paper, we have examined and compared no less than 20 
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Figure 1. NGC 891 from the Digitized Sky Survey. The num- 
bered red boxes indicate the locations in which the flux was mea- 
sured to create the spectra. 

independent data sets 2 ranging in wavelength from A 1.2 /jm 
to 850 /im using both space-based and ground-based instru- 
ments (Sect. 2). To the limits of the variety of spatial res- 
olution, coverage, and sensitivity achieved, we compare the 
extraplanar extent of the emission and, most importantly, 
we derive the SED at different disk and extraplanar loca- 
tions. 

Because of the proximity of NGC 891, its almost ex- 
actly edge-on orientation, its current active star formation 
(SF) and its prominent extraplanar emission (see Sect. 1.1), 
this galaxy has been a prime target for studies of the disk- 
halo interaction. Although the specific vertical height, z, at 
which gas might be considered in the halo is debatable, a 
conservative estimate is 1 kpc (21.5 arcsec). Since all of the 
20 data sets examined have higher resolution than 21.5 arc- 
sec (see Sect. 2), it is possible to distinguish the disk from 
extraplanar (halo) region in the observations presented in 
this paper. 

In the next section we describe NGC 891. In Sect. 2, we 
discuss the observations and data reduction. The morphol- 
ogy of NGC 891 is discussed in Sect. 3, and the SED results 
are given in Sect. 4. Sects. 5 and 6 present the discussion 
and conclusion, respectively. 

1.1 NGC 891 

NGC 891 (Figure 1, Table 1) is believed to be similar to the 
Milky Way (MW) in B-band luminosity, rotational velocity, 
and Hubble type. It is 2.5 times brighter in the IR, and may 
have twice as much gas and dust than the MW (see Scoville 
et al. 1993), although the difference is not strong, given un- 
certainties. There has been much discussion in the literature 
as to whether or not NGC 891 is a 'starburst' galaxy, since 

2 A A 160 image was also examined, but found to be of too 
low spatial resolution to be useful. 
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Table 1. Basic Galaxy Parameters of NGC 891 



parameter 


value 


reference 


Hubble type 


SA(s)b 


RC3 a 


RA(J2000) 


02h 22m 33.4s 


NED b 


DECf 12000 "I 


+42° 20' 57" 


NED 


rcdshift 


528 km/s 


NED 


distance 


9.6 Mpc 


Strickland et al. (2004) 


optical major axis 


13.5 arcmin 


RC3 


optical minor axis 


2.5 arcmin 


RC3 


major axis angle 


22° 


NED 


inclination angle 


89° 


Baldwin & Pooley (1973) 


rotational velocity 


225 km/s 


Sofuc (1997) 


B band luminosity 


7.8 x 10 9 L© 


RC3 


SFR 


3.8 M Q /yr 


Popescu et al. (2004) 


M gaa (H I & H 2 ) 


4.6 X 10 9 M 


Guelin et al. (1993) 






Dupac et al. (2003) 


Mdust 


1.9 - 7 x 10 7 M C 


Alton et al. (2000) 






Popescu et al. (2004) 






Galliano et al. (2008a) 


MpAH 


3.3 x 10 6 M Q 


Galliano et al. (2008a) 



a RC3 = Third Reference Catalogue by de Vaucouleurs et al. (1991) 
b NED = NASA Extragalactic Database 
(http://nedwww.ipac . caltech.edu/) 

c The range reflects different results from various authors. 

its star formation rate (SFR) appears to be intermediate be- 
tween that of a normal spiral and a starburst galaxy; it may 
be in a phase that saw more powerful starbursting in the 
past (see Temple et al. 2005). NGC 891's current SF activ- 
ity is clearly indicated by a variety of measures, including the 
detection, in Ha, of many bubbles, shells, and supershells in 
the mid-plane of this galaxy (Rossa et al. 2004). 

The molecular g cLS, ELS observed in CO, resembles that 
of the MW and is mainly confined to a thin disk (400 pc 
thickness) (Scoville et al. 1993) although a halo component is 
also observed (see below). Emission from cold dust (average 
temperature of ~ 18 - 24 K; Dupac et al. 2003), as observed 
at A 1.3 mm, 850 /im and 450 /im, correlates spatially with 
the molecular gas in the disk (Alton et al. 1998; Israel 1999). 
Spatial correlations have also been found between cold dust 
at A 850 /im, the 7.7 /im PAH feature and very small grain 
(VSG) emission observed at A 14.3 /im that is not specifically 
associated with starburst regions (Haas et al. 2002). 

As indicated in Sect. 1, NGC 891 has extended extra- 
planar gas in a variety of ISM components. The galaxy has 
widespread diffuse ionized gas (DIG) in its halo (Rossa & 
Dettmar, 2003; Rossa et al. 2004) up to 2.2 kpc above the 
galactic plane. It is one of the brightest of all nearby edge-on 
spirals, suggesting pronounced disk-halo mass transfer (Al- 
ton et al. 2000). The Ha halo emission is more prominent 
and extended in the z direction on the northern side of the 
galaxy than the southern. Some of this asymmetry may be 
explained by dust obscuration (Kamphuis et al. 2007), but 
Rossa et al. (2004) interpret this as being due to a higher 
SFR in the northern part of the disk than in the southern 
part. The assymetry is also seen in the radio continuum ob- 
servations (Dahlem et al. 1994) and these are not affected 
by dust. They therefore, support the latter interpretation. 

The broadest halo component is that of H I, which also 
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shows a north-south asymmetry in disk thickness (Swaters 
et al. 1997), and is now known to extend up to at least 22 kpc 
above the galactic plane (Oosterloo et al. 2007). Also, the H 

1 halo is asymmetric: It is clearly thicker in the north than in 
the south of the galaxy (Swaters et al. 1997). Molecular gas, 
as observed in CO, has been detected to 1 - 1.4 kpc above the 
plane (Garcia-Burillo et al. 1992). As for dust in the halo, 
optical absorption features have been detected to heights of 

2 kpc (Howk & Savage 2000). Burgdorf et al. (2007) have 
also detected A 16 and A 22 fira emission extending up to 5 
kpc from the midplane. Recently, Rand et al. (2008) have 
obtained IR spectra showing the [Nell] and [Nelll] MIR lines 
of two positions in NGC 891 's halo which suggest that the 
radiation field may harden with height from mid-plane. 



2 OBSERVATIONS AND DATA REDUCTION 

This paper makes use of both archival data that we have 
reduced as well as previously published data. In this section, 
we first describe the data that we reduced ourselves and then 
the origin of the data sets that were obtained from others. 

2.1 ISO Observations and Data Reduction 

The Infrared Satellite Observatory (ISO) archive was ac- 
cessed to obtain two previously unpublished data sets that 
were acquired in 1997 and 1998. The data were taken with 
ISO's camera, ISOCAM, using the long wave channel, with 
narrow band filters, LW1, LW4, LW5, LW6, LW7, LW8, 
LW9, (centered at A 4.5, 6.0, 6.8, 7.7, 9.6, 11.3, 14.9 £tm, re- 
spectively), and wide band filters, LW2 and LW3 (A 6.7 and 
14.3 lira, respectively). All of these observations were with 
the CAM01 observing mode, which is the standard mode for 
photometric imaging (Blommaert et al. 2003). 

The galaxy was not completely contained in the field of 
view. Only 3.2 arcmin x 3.2 arcmin (8.9 x 8.9 kpc) were 
viewed at a time, and the galaxy was not always centered in 
the field. There were multiple pointings which had to be mo- 
saiced together (see below) to form the maps shown in Sect. 
3.1. The raw data consist of a sequence of 2-dimcnsional 
frames taken at successive times for any given wavelength 
band, forming a 'cube' of data. Table 2 summarizes the ob- 
servation parameters. 



Table 2. ISO observing and map parameters. 



Parameter 


LW1 


LW4 


LW2(wide) 


LW5 


wave bands 
LW6 


LW7 


LW8 


LW3(wide) 


LW9 


central wavelength a (/im) 


4.5 


6.0 


6.7 


6.8 


7.7 


9.6 


11.3 


14.3 


14.9 


wavelength range" (/im) 


4.0-5.0 


5.5-6.5 


5.0-8.5 


6.5-7.0 


7.0-8.5 


8.5-10.7 


10.7-12.0 


12.0-18.0 


14.0-16.0 


TDT number'' 


65101201 


65101201 


65101201 


65101201 


65101201 


65101201 


65101201 


65101201 


65101201 








84101766 










84101766 




pixel field of view*" (arcsec) 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


PSF^ fFWHlVH (arapc) 


7.3 


7.3 


7.3 


7.9 


7.9 


7.9 


7.9 


7.9 


8.4 


date of observation 


28 Aug 97 


28 Aug 97 


28 Aug 97 


28 Aug 97 


28 Aug 97 


28 Aug 97 


28 Aug 97 


28 Aug 97 


28 Aug 97 








5 March 98 










5 March Q8 




no. of on-source pointings 6 


7 


7 


15 


7 


7 


7 


7 


15 


7 


mean # frames/pointing 


16 


16 


13 


16 


13 


16 


16 


16 


16 


integration time per frame (s) 


10.08 


10.08 


10.08 


10.08 


10.08 


10.08 


10.08 


10.08 


10.08 


total on source time^ (min) 


18.8 


18.8 


36.8 


18.8 


15.3 


18.8 


18.8 


40.3 


18.8 


sky coverage (arcmin) 


3.2 


3.2 


3.2 


3.2 


3.2 


3.2 


3.2 


3.2 


3.2 


RMS of final maps (MJy/sr) 


0.12 


0.14 


0.13 


0.22 


0.16 


0.10 


0.28 


0.19 


0.19 


dynamic range of final maps^ 


189 


202 


370 


167 


600 


119 


138 


136 


91 



a Blommaert et al. (2003) 

^ a unique number that identifies the observation. 
c See text Sect. 2.1 

d The point spread function (PSF) full width at half maximum (FWHM). 
e This is the number of different pointing positions on-sourcc. 

f Total integration time per on-sourcc pointing times the number of on-source pointings. 
9 maximum of map divided by the RMS (noise) of final map. 
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For the data reduction, we used the Cam Interactive 
Reduction (CIR) package (Chanial 2003) which runs in In- 
teractive Data Language (IDL). A thorough review of the 
data reduction procedures can be found in Galliano (2004) . 

We first performed a dark correction following the char- 
acterization of Biviano et al. (1998). A second-order correc- 
tion is also carried out during this step, depending on the 
detector temperature and the time since activation, and a 
short-drift correction is applied. 

Cosmic rays were then removed by running an auto- 
matic multi-resolution de-glitching routine which masks the 
glitches that are short in duration and have an intensity 
above 6u (where a is the noise level). Further manual de- 
glitching was then carried out, frame by frame, to remove 
the remaining glitches that did not meet these criteria. These 
were typically pixels that were immediately adjacent to ones 
that had been automatically deglitched. 

Memory effects caused by non-instantaneous temper- 
ature changes in the photoconductors were then corrected 
using the Fouks-Schubert (1995) method and the Wozniak 
implementation (see Galliano 2004). 

We removed a narrow border of each frame, since the 
pixels at the very outer edges of the CCD have higher noise. 
This step has not significantly reduced the field of view be- 
cause the various individual frames overlap spatially and the 
final field of view is dictated by the mosaic of these frames. 

The data were separated into blocks of frames corre- 
sponding to each wave band. The blocks were then regis- 
tered to the same sky position and medianed together to 
create a single image at each wave band. 

We then flat fielded each image using flat fields main- 
tained in a flat library at the ISO Data Center and also 
corrected for flat distortion which is a known ISOCAM er- 
ror. At this point, the data were calibrated into mjy/pixel 
units using the calibration terms given by Blommaert et al. 
(2003). 

Background emission was measured in a region that was 
far from the main emission (greater than 15 kpc). This back- 
ground emission was then subtracted from the image. 

Finally, we created error maps for each wave band (not 
shown), which represent pixel- by-pixel RMS values derived 
during cube contraction, propagating through errors in the 
various processing steps. 

The remaining source of error not included in the error 
maps is the ISOCAM absolute flux calibration error, which 
does not affect the pixel-by-pixel variation within a given 
image. The calibration error is by far the largest of the er- 
rors, at about 20% (Pagani et al. 2003; Coia et al. 2005). 

The effect of a dead pixel column on the detector re- 
mains on all of our ISOCAM images. This region is blanked 
in Figures 2, Al, A3 and A4. We do not make any measure- 
ments in regions affected by this known error. 

2.2 Spitzer Observations and Data Reduction 

We use Spitzer Space Telescope data with both the Multi- 
band Imaging Photometer (MIPS) as well as the Infrared 
Array Camera (IRAC). Observing and map parameters for 
the Spitzer data can be found in Table 3. 

MIPS scan map data were only initially available for As- 
tronomical Observation Requests 14815744, 14815488, and 
14816000, which were data taken as part of an observing 



program led by C. Martin, so only those data were used in 
our analysis. The raw data were processed and combined 
into final mosaics using the MIPS Data Analysis Tools ver- 
sion 3.06 (Gordon et al. 2005). The data processing for the 
A 24 /im data differs significantly from the A 70 /im data, so 
the data reduction procedures are given separately. 

The individual A 24 pirn frames were first processed 
through a droop correction (removing an excess signal in 
each pixel that is proportional to the signal in the entire 
array) and were corrected for non-linearity in the ramps. 
The dark current was then subtracted. Next, scan-mirror- 
position dependent flats were created from the data in each 
Astronomical Observation Request (AOR) and were applied 
to the data to remove dark spots in the image that are caused 
by contamination on the MIPS scan mirror. Scan-mirror- 
position independent flats were created from the data in each 
AOR and were applied to the data to remove negative latent 
images from observations of bright objects that preceded the 
individual AORs and to correct sensitivity variations across 
the detector. Following this, the zodiacal light emission was 
fit with a third-order polynomial function in each scan leg, 
and the function was then subtracted from the data. A ro- 
bust statistical analysis was then applied in which the values 
of cospatial pixels from different frames were compared to 
each other, and statistical outliers such as cosmic rays were 
masked out. After this, a final mosaic was made with pixel 
sizes of 1.5 arcsec, any residual background in the image was 
subtracted, and the data were calibrated into astronomical 
units using the calibration term 0.0447 MJy sr^ 1 MIPS~„ it . 
The uncertainty in the calibration is 4%. While the FWHM 
of the PSF is ~ 6 arcsec, significant extended emission from 
the Airy rings can be found outside the central peak. 

In the A 70 /im data processing, the first step was to fit 
ramps to the reads to derive slopes. In this step, readout 
jumps and cosmic ray hits were also removed, and an elec- 
tronic nonlinearity correction was applied. Next, the stim 
flash frames taken by the instrument were used as respon- 
sivity corrections. The dark current was subtracted from the 
data, and an illumination correction was applied. Short term 
variations in the signal (often referred to as drift) were re- 
moved from the data, which also subtracted the background 
from the data. Next, a robust statistical analysis was applied 
to cospatial pixels from different frames in which statistical 
outliers were masked out. Once this was done, final mosaics 
were made using square pixels of 4.5 arcsec. The residual 
backgrounds were measured in regions outside the optical 
disk of the galaxy and subtracted, and then the flux cal- 
ibration factor of 0.2 MJy sr" 1 MIPS^ it was applied to 
the data. The uncertainty in the calibration term was 7%. 
The final image is very strongly affected by bright and dark 
streaks outside the plane of the galaxy that are related to 
latent image effects. Also, while the FWHM of the PSF is 
~ 18 arcsec, significant extended emission from the Airy 
rings can be found outside the central peak. 

Aside from uncertainties in the MIPS calibration terms 
indicated above, there are additional non-linearities and 
other effects which we estimate puts errors on flux mea- 
surements, when comparing with other bands, of order 10 
to 15%. 

The IRAC observations at A 3.6, 4.5, 5.8, and 8.0 /im 
were taken from the Spitzer archive. These are the post- 
basic calibrated data (data that are already reduced and cal- 
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ibrated) in four MIR wave bands: A 3.6, 4.5, 5.8, and 8.0 /Ltm. 
These data were part of the Brown Dwarf Galaxy Haloes 
observing program led by Giovanni Fazio (Program ID #3). 
Calibration uncertainties are of the order of 5% (Reach et 
al. 2005). 

Additional data reduction on these IRAC data were 
as follows: First, a sky subtraction was performed on each 
band. Then, maximum 'infinite aperture' corrections of 0.94, 
0.94, 0.78, and 0.74 have been recommended for the A 3.6, 
A 4.5, A 5.8, and A 8.0 wavebands, respectively (Reach et al. 
2005) for extended sources. Since these corrections are of or- 
der of the absolute calibration error for the 3.6 and 4.5 mi- 
cron bands, we make no corrections for those bands. For the 
longer wavelength bands, however, we have applied a correc- 
tion for all disk emission flux measurements in the amount 
of 0.94 and 0.78 at A 5.8 and A 8.0, respectively. These val- 
ues assume that the maximum error must be applied in a 
direction along the disk whereas the disk is well defined and 
finite in the z direction (see Figure 9 plus the Spitzer web- 
site ). Since infinite aperture corrections are unknown for 
the halo regions, we omit the 5.8 and 8.0 /im data points in 
the SED halo calculations (Section 4.3). The displayed maps 
are shown without this correction. 

For additional IRAC details and map rms, we refer to 
Table 3. 

2.3 2MASS and JCMT Observations 

The remaining data sets include the 2 Micron All Sky Sur- 
vey (2MASS) images in the J, H and K bands, which were 
obtained from the Skyview archive 4 (Jarrett et al. 2003), 
as well as the Sub-millimeter Common User Bolometer Ar- 
ray (SCUBA) data of the James Clerk Maxwell Telescope 
(JCMT) at A 450 and 850 /im which were kindly supplied by 
E. Xilouris (see Alton et al. 1998). Calibration uncertainties 
are approximately 2% for 2MASS (Cohen et al. 2003) and 
15 - 25% for SCUBA (Alton et al. 1998 and Tilanus 2005). 

For the SCUBA data, we adopt the estimates of Israel 
et al. (1999) who indicate that about 5% CO(3-2) contam- 
ination may be in the A 850 /tm band with < 10% toward 
the nuclear region, and that negligible CO (6-5) contamina- 
tion is in the A 450 /xm band. Since these are less than the 
calibration uncertainties, we do not correct for CO in these 
bands. 

For observation and map parameters for these data, see 
Table 4. 

2.4 Final Maps 

For comparative purposes, all data were regridded to the 
same pixel size (1 arcsec) and the flux converted to units that 
could be compared, map to map. The units are therefore, 
either mjy arcsec -2 or MJy sr -1 (note that the conversion 
from mjy arcsec -2 to MJy sr -1 requires a multiplication by 
42.55). Most data in this paper are presented in MJy sr -1 
units, though we have also retained the original Jy beam -1 
units for the JCMT/SCUBA data to present those maps in 
a form that is standard at those wavelengths. 

3 http : //ssc . spitzer . caltech . edu/irac/calib/ extcal/ 

4 http : //skyview . gsf c . nasa . gov/cgi-bin/ skvadvanced . pi 



Table 3. Spitzcr observing and map parameters. 



Parameter 


C ' i /T"D A f 1 

bpitzcr / 1KAG 








bpitzcr/MLFb 




Tpfpvpnpp wfivpIpti or Vi ^ (imn} 

1 i- 1 1 1 WOiVLUvll^Lrll \ Lit 111 / 


3.6 


4.5 


5.8 


8.0 


24 


70 


central wavelength (lira) 


3.6 


4.5 


5.8 


8.0 


23.68 


71.42 


bandwidth^ (pm) 


1.0 


1.0 


1.5 


3 


8 


20 


pixel field of view (arcsec) 


1.2 


1.2 


1.2 


1.2 


1.5 


4.5 


PSF (FWHM) (arcsec) 


1.9 


2.0 


1.9 


2.2 


6 


18 


date of observation 


8 Sept 2004 


8 Sept 2004 


8 Sept 2004 


8 Sept 2004 


4 Sept 2005 


4 Sept 2005 


RMS of final maps (MJy/sr) 


0.026 


0.018 


0.023 


0.057 


0.048 


0.050 


dynamic range of final maps 


136 


169 


309 


611 


587 


490 



a See Table 2 for explanation of parameters that are not described here. 

^ The reference wavelength is the name for the wave band that is often used in the literature. 
c The central wavelength is the actual center of the wave band. 
^ Wavelength range of filter. 



Table 4. 2MASS and JCMT/SCUBA observing and map parameters. 



Parameter 


2MASS 






JCMT/SCUBA 




waveband (/an) 


J 


H 


K 


450 fj,m 


850 iim 


central wavelength (fJ.ro) 


1.24 


1.66 


2.16 


442 


862 


bandwidth^ (pm) 


0.162 


0.251 


0.262 


50 


130 


PSF (FWHM) (arcsec) 


1 


1 


1 


7.5 


15.7 


RMS of final maps 


1.7 X 10~ 7 MJy/sr 


3.1 X 10~ 7 MJy/sr 


2.8 X 10~ 7 MJy/sr 


0.83 Jy/beam 


0.010 Jy/beam 



f 

as 
3 

a Name of wave band. »>. 

a. 
Co 
s 

00 
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3 RESULTS 

3.1 IR Images and Emission from the Disk 

Figures 2-4 and Al - A4 show the MIR surface brightness 
of NGC 891 in both contours and grey scale. The ISO A 7.7 
/xm narrow wave band and the IRAC A 8 /im wave band are 
particularly good at selecting PAH spectral features and are 
shown in Figure 2 (with the other ISO PAH bands shown 
in Figure Al). The A 4.5 /jm wave band is shown in Figure 
3, and the A 5.8 /im wave band is shown in Figure 4. Figure 
A4 shows the two wide bands from ISO. In each figure the 
foreground stars have been indicated by small crosses 5 and 
the dead pixel column is contained within the blanked white 
region of the ISO images. 

Figures 2-4 (and Al - A4) show that most of the emis- 
sion is coming from the disk of the galaxy, especially from 
the location of the central dust lane. These surface bright- 
ness maps show prominent peaks toward the nucleus and 
secondary maxima on either side of the center. These sec- 
ondary peaks are associated with spiral arms and bright H 
II regions (Kamphuis et al. 2007). Figure 3 shows that the 
contours at A 4.5 fim form a bulge at the center, whereas at 
longer wavelengths, we do not see this bulge. Comparison of 
major axis profiles (not shown) also reveal substantial dif- 
ferences between the short and long wavelength emission. 
This change in morphology with increasing wavelength il- 
lustrates the changing dominance of spectral components. 
The A 3.6 /jm wave band contains the 3.3 fim PAH feature 
and Le Coupanec et al. (1999) concluded that the A 4.5 /jm 
flux in NGC 891 contained a strong warm dust component. 
However, Figure 5, which shows the emission contours in 
these two wave bands overlaid on a stellar map (2MASS K 
band), shows an excellent correspondence between all three 
bands. This suggests that the A3.6/im and 4.5 /im bands 
are dominated by starlight, a result that is borne out by our 
SED results (Sect. 4.3). We are therefore in agreement with 
most of the results from SINGS and other Spitzer studies, 
which have demonstrated that the 3.6 and 4.5 /j,m emission 
is almost entirely stellar in origin (Draine et al. 2007b; Smith 
et al. 2007). 

Figures 6 and A5 show the warm to cold dust emission 
as wavelength increases from A 24, through 70, 450, and 850 
/im. Comparing these maps with Figures 2-4 and Al - 
A4, we find that, along the plane of NGC 891, the local 
maxima of emission for PAH-sensitive bands are at simi- 
lar locations as the local maxima for both warm dust and 
cold dust. Therefore, to first order, PAHs and dust emission 
correlate in the disk. This result is likely simply a state- 
ment that where there is a density enhancement of dust, 
there will also be a density enhancement of molecular mate- 
rial and PAHs as well. The dust temperature and source of 
PAH excitation, however, impose additional important con- 
straints on the apparent distribution of these components. 
PAHs and cold dust are known to be correlated from previ- 
ous work (Haas et al. 2002; Bendo et al. 2008), suggesting 
that PAH excitation is possible from a widespread stellar 
FUV-emitting population, rather than in SFRs alone. 

As for warm and cold dust, we investigate the relation- 
ship between these components via the A 24 /im to 850 /jm 
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Figure 2. Surface brightness contours of the ISO A 7.7 (im and 
the IRAC A 8 /im wave bands that select PAH emission. The 
crosses denote foreground stars. The contours are 0.4 (3<r), 1, 3, 
5, 8 MJy/sr for A 7.7 /tin, and 0.17 (3a), 0.4, 1.3, 5, and 7 MJy/sr 
for A 8 (im. 



5 Foreground stars were identified using the 2MASS K band map. 
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Figure 3. Surface brightness contours of the IRAC A 4.5 (im 
wave band that selects mostly the stellar continuum. Symbols as 
in Figure 2. Contours arc 0.05 (3a), 0.15, 0.4, 1.2, 2.5, and 5 
MJy/sr. 

surface brightness ratio map, after first smoothing the A 24 
/xm map to the same resolution as the A 850 /im map. The 
result is shown in Figure 7. From this figure, it seems that 
the cold dust emission (A 850 /im) is slightly more extended, 
since the ratio tends to decrease outward. This agrees with 
our results for the z distribution as well (see Sect. 3.3). In 
addition, in the disk there are two "blobs", or maxima, in 
the ratio, on either side of the centre identifying regions in 
which the dust is hotter than average. The southern maxi- 
mum is located at the same position as the surface brightness 
maxima seen in the previous maps (Figures 2-5 and A5). 
The northern ratio maximum is located further away from 
the centre, at a declination of 42° 22' 45". Both these max- 
ima are located in regions in which the Ha emission is also 
prominent (e.g. see the Ha map in Kamphuis et al. 2007). 
The dust heating is likely related to in-disk star forming 
regions at these locations. 

3.2 Major Axis Profile 

On large scales, Swaters et al. (1997) found a north-south 
asymmetry in the HI emission, such that the southern emis- 
sion is more extended radially along the major axis. Popescu 
& Tuffs (2003) subsequently found a similar result in the 
A 200 emission. Asymmetries have also been seen in 
the molecular distribution (see Sofue & Nakai 1993). Given 
the correlations between bands indicated above, we inves- 
tigate whether any global asymmetries are present in our 
data using a representative PAH band which we take to be 



Figure 4. Surface brightness contours of the A 5.8 (im wave band 
that selects the MIR continuum. Symbols as in Figure 2. Contours 
arc 0.068, 0.17, 0.6, 2, 3.4 MJy/sr. 

the Spitzer IRAC A 8 /im data (Figure 2), and the Spitzer 
MIPS A 24 /im data (Figure 6) since the entire galaxy has 
been imaged in these wavebands (as opposed to the narrower 
band but truncated ISO images of Figure 2). The result is 
shown in Figure 8. Since the noise is very low (see figure cap- 
tion), virtually all the structure seen in this profile is real. 
The broadscale distribution of emission is consistent with 
what has been seen at other MIR wavelengths (e.g. Dupac 
ct al. 2003) and this PAH profile also compares well with the 
lower spatial resolution ISOPHOT observations (Mattila et 
al. 1999, their Figure 5a). However, to our knowledge, Fig- 
ure 8 presents the highest resolution and highest S/N major 
axis profile yet obtained. 

The profile shows a peak at the centre and two broad 
peaks on either side which, together form a broad 200 arcsec 
central plateau of PAH emission. This emission contains the 
most active part of the disk encompassing the nucleus and 
the secondary peaks at -40 arcsec (south of the nucleus) 
and +40 arcsec to +60 arcsec (north of the nucleus). Much 
substructure is also visible. As for asymmetries, we can say 
that there is some enhancement on the northern side within 
a radius of w 150 arcsec and a minor enhancement towards 
the south at radii larger than this. 

3.3 The Vertical Distribution & Extraplanar 
Emission 

The 'edges' of the emission shown in Figures 2-4 and Al 
- A4 are not smooth and show varying widths as a function 
of location in the disk. Discrete extensions, such as seen in 
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Figure 5. IRAC A 3.6 and 4.5 (im emission contours overlaycd 
on the K band image. The contours are the same as for Figure 3. 

other galaxies (e.g. Irwin et al. 2007) or in the DIG of this 
galaxy (Rossa & Dettmar 2003, see Sect. 1.1) are not promi- 
nent at the S/N of these data, although there are hints of 
such discrete features at slightly lower S/N. We know that 
vertical discrete features are indeed present in the dust dis- 
tribution from the optical absorption observations of Howk 
& Savage (1997) and the Ha observations of Rossa et al. 
(2004). 

To further investigate the presence of extraplanar dust 
in these data, we have averaged the emission parallel to the 
major axis (within ± 2.5 arcmin radially from the minor 
axis) to create an average vertical profile. While such av- 
eraging removes information about any discrete features, it 
does increase the S/N, improving the possibility of detecting 
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Figure 6. MIPS A 24 (im and 70 (im maps. The contours start 
at 3cr. Symbols as in Figure 2. Contours are 0.145, 0.3, 1, 3, 10, 
50 MJy/sr for A 24 /im, and 1.5, 3, 7, 15, 40, 120, 350 MJy/sr for 
A 70 /tin. 
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Figure 7. The A 24 (im to 850 fim surface brightness ratio map. 
The contours are 2, 3, 4, 5, 6, 8, 9.5. 
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Figure 8. Surface brightness profiles along NGC 891's disk for 
the A 8 fim band (rms noise of 0.025 MJy sr — x ) and the A 24 
fim band, obtained from an average of ± 150 arcsec in z above 
and below the plane of the image in Fig. 2. Positive along the x 
axis corresponds to the northern side of the disk and negative the 
southern side. The A 24 fim flux has been normalized to the A 8 
fim data and can be can be recovered by multiplying by 333. The 
A 24 fim noise, prior to normalization, is 0.86 MJy sr" 1 . 



widespread high latitude dust. We have also taken a subset 
of the data sets and smoothed them to a common spatial 
resolution (either 7.9 arcsec or 20 arcsec) so that we can 
compare the vertical extent of the emission in various bands. 
Results are shown in Figure 9, in which a typical profile (in 
this case the narrow band A 7.7 /im ISO PAH band) is shown 
in linear scale to emphasize the exponential fit, and in a log 
scale in Figures 10 and 11 which compare the various pro- 
files at 7.9 arcsec and at 20 arcsec resolution, respectively. 



As was the case for the major axis profile, these average z 
slices have very small vertical error bars (see Figure cap- 
tions), and have been omitted from the figures for clarity. 
The main error, for the purpose of comparing the width of 
the emission in the various wavebands, is in the horizontal 
direction, associated with the spatial resolution (±4 arcsec 
for Figures 9 and 10, and ± 10 arcsec for Figure 11). 

Figure 9 is a typical MIR vertical profile and shows 
that an exponential fall-off of emission fits the data well in 
both the core and wing region. This is also true of the other 
MIR bands (fits not shown, but see also the exponential 
fits of Burgdorf et al. 2007 and Rand et al. 2008). We have 
illustrated the narrow band A 7.7 fim ISO profile in partic- 
ular, since it is excellent at isolating PAH emission. Note 
that the wider waveband A 8 fim IRAC profile matches the 
A7.7 fim ISO profile to within errors (cf. Figure 10). Single- 
component exponentials fit all data sets well, although some 
departures are seen at very low levels. The exponential scale 
heights that result from these fits, z e , are not sufficiently 
different from the spatial resolution to be convincingly pre- 
sented. However, of the dust-related data sets at their orig- 
inal un-smoothed resolution, the IRAC A 5.8 fim and IRAC 
A 8.0 fim data sets have high enough spatial resolution (Ta- 
ble 3) for a z e measurement. For these, we find (over the 
same 5 arcmin averaged region as Figure 9) z e = 302 ± 70 
pc and z e = 248 ± 50 pc at A 5.8 and 8.0 fim, respectively, 
where the error bars are dominated by the variation between 
the north-west and south-east sides of the galaxy. Our scale 
height for the PAH-dominated A 8.0 fim band over this av- 
eraged region is slightly lower than the single-component 
z e estimate of 330 - 530 pc found for A 17 fim PAH emis- 
sion measured at a single position by Rand et al. (2008). 
Given the possible effects of unknown in-disk extinction on 
the vertical scale height (see Rand et al. 2008) as well as 
uncertainties related to aperture corrections for the IRAC 
data (Sect. 2.2), two-component fits are not warranted. 

Figure 10 shows a comparison between bands at high 
spatial resolution, for which the error in z is ± 4 arcsec, or 
0.2 kpc. The logarithmic scale emphasizes lower levels for 
which there are some departures from perfect exponentials. 
The broadest component is clearly in K band, due to the 
dominance of starlight in this band and the presence of the 
galaxy's bulge. The cool dust component, illustrated by the 
A 450 fim emission, is the widest dust band and is measured 
to z m 29 arcsec (1.3 kpc) at 3<r (average z of both sides 
of the vertical profile). The PAH emission (A 7.7 and 8 fim) 
and warm dust emission (A 24 fim) have z widths that agree 
within errors. Also, their measurable maximum extents are 
to z = 49 arcsec (2.3 kpc), 53 arcsec (2.5 kpc), and 54 arcsec 
(2.5 kpc) at A 7.7, A 8 fim, and A 24 fim, respectively, again 
agreeing within errors. Although the maximum extent of the 
various components depends on the varying S/N, it is clear 
that each of these components reveal emission at z in excess 
of 1 kpc. It is therefore clear that NGC 891 has extraplanar 
emission in warm dust, cool dust, and PAHs. Kamphuis et 
al. (2007) have also recently examined the MIPS A 24 fim 
emission from NGC 891 and find emission to z of 2.3 kpc, 
in agreement with our results. Burgdorf et al. (2007) have 
also found extraplanar dust emission extending to z =5 
kpc, at a specific location rather than an average, and at 
the wavelengths, A 16 and A 22 fim. 

As for the lower resolution profiles (Figure 11), the spa- 
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Figure 9. Surface brightness profile (left) for the A 7.7 fim wave band (red) as a function of z height, with an exponential fit (green). 
We have zoomed into the wing region (right) to better see the fit. The spatial resolution of the data is 7.9 arcsec. The error bars have 
been omitted for clarity (see text). The vertical errors are ± 0.028 MJy sr — 1 . 
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Figure 10. Surface brightness profiles for the A 7.7 (im (red, as 
in Figure 9), 2MASS K band (green), A 24 fim MIPS band (blue), 
A 450 fim SCUBA band (magenta), and A 8 fim IRAC band (cyan) 
as a function of z height; all data have been smoothed to 7.9 arcsec 
resolution prior to forming the profile. All data have been cut off 
at the 3cr level which, prior to scaling, is 0.085, 0.033, 1.6, 2.6, 
and 0.13 MJy sr _1 at A7.7 /im, K band, A24 /mi, A450 fim, and 
A 8 fim, respectively. The horizontal error bars are ± 4 arcsec. 
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Figure 11. Surface brightness profiles for A 7.7 /an (red), H I 
emission from Swaters (1997) (green) and A 850 (im SCUBA band 
(blue) as a function of z height. These data have been smoothed to 
20 arcsec resolution and have been cut off at the 3<r level (with the 
exception of the HI). The horizontal error bars are ± 10 arcsec. 



tial resolution (20 arcsec = 0.9 kpc) is less useful for exam- 
ining extraplanar gas quantitatively. However, it is apparent 
that the H I distribution is widest (see also Sect. 1), followed 
by cold dust at A 850 fim and PAHs. 

4 SPECTRAL ENERGY DISTRIBUTION 
4.1 Constructing the SED 

We have constructed spectra using the flux from different re- 
gions of the galaxy in each of the wave bands. The adopted 
boxed regions are shown in Figure 1 and described in Ta- 
ble 5. They represent a potential diversity of environments 
within which the S/N of the data was sufficient at all wave- 
bands for good results. Each box is 18 arcsec x 18 arcsec 
(838 pc x 838 pc) in size, a value chosen because all images 



were smoothed to 18 arcsec resolution to match the low- 
est resolution A 70 fim data before the measurements were 
made. The four in-plane boxes represent the centre of the 
galaxy and three disk regions. The centre contains the disk 
that passes along our line of sight (old and young stars and 
ISM), old stars from the bulge and the nucleus. The disk 
region in the southern half is located at a local maximum. 
Although there may be some contributions from the bulge to 
the northern and southern disk regions, our results for these 
regions should be strongly dominated by the disk compo- 
nents that are along the line of sight. All in-disk regions fall 
within the broad plateau of emission shown in Figure 8 in- 
cluding the outer north disk at a radial distance of 96 arcsec 
(4.47 kpc). Finally there are two halo regions, each centered 
between 1 and 2 kpc above the plane. These regions should 
therefore provide one of the first MIR spectra of halo ISM 
in an external galaxy (with Rand et al. 2008 being the only 
other one at the time of this publication). 
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Table 5. Regions of NGC 891 Chosen for Spectral Analysis 





name of 


RA & DEC a 


offset (R, zf 


# 


region 


(h m s), (° ' ") 


(kpc, kpc) 


1 


centre 


2 22 33.092, 42 20 50.26 


0, 


2 


south disk 


2 22 30.928, 42 20 05.26 


2.51, 


3 


inner north disk 


2 22 35.258, 42 21 52.26 


2.79, 


4 


outer north disk 


2 22 37.153, 42 22 30.26 


4.47, 


5 


central halo 


2 22 30.927, 42 21 06.26 


0, 1.40 


6 


halo over disk 


2 22 33.002, 42 21 59.26 


2.79, 1.26 



0.046 



M d 



(1) 



and the ratio of ionized PAHs to total PAHs (fp A H+ ) 
M PAH + 



fpAH+ 



MpAH° + M PAH + 



(2) 



a Coordinates refer to J2000.0 

k Distance from the center of the galaxy along and perpendicular 
to the major axis, respectively. 

For the flux measurements, the error bars are due to the 
published flux errors given for each instrument. To summa- 
rize from Sect. 2, these are approximately 20% for ISOCAM, 
5% for IRAC, 4% - 15% for MIPS, 2% for 2MASS, and 15 - 
25% for SCUBA. These are the errors that dominate when 
comparing between different wave bands. We present the 
spectra fitted by SED models in Figure 12. Note that we 
only use upper limits for the 70 micron measurements of 
the halo because the extraplanar region is severely affected 
by latent image effects that create streaking in the images. 
This data point is carrying lower weight in the SED mod- 
elling. The next section describes this model. 



4.2 Modeling the SED 

The model we use fits the observed SEDs for a situation 
in which a variety of dust components with realistic optical 
properties are exposed to an incident radiation field from 
a mixture of stars, all within the same beam. The model 
includes the following IR emission sources: 

1) PAHs - which can be singly ionized (PAH + ) or neutral 
(PAH°), 

2) Dust - which include VSGs and large grains composed of 
graphites and silicates, and 

3) Old Stars 

We adopt the dust size distributions derived by Zubko 
et al. (2004) by fitting the diffuse Galactic ISM emission, 
extinction and depletion patterns, with PAH, graphite and 
silicate grains, and assuming solar abundances. Thus we are 
assuming that the instrinsic dust properties in NGC 891 are 
similar to the ones in the MW. In fact, this is a good as- 
sumption because both galaxies are solar metallicity spiral 
galaxies (12+log(0/H)=8.9, Otte et al. 2001). The PAHs are 
assumed to have the optical properties (e.g. absorption cross 
sections) from Draine & Li (2007a). Note that, for the pur- 
pose of the model, PAHs will be included in the 'dust' com- 
ponent, rather than considered as 'large molecules', given 
the similarity between the largest PAHs and smallest car- 
bon grains. Thus the model parameter, M d , represents the 
total dust mass; this includes PAHs, graphites and silicates. 

The PAH component introduces two parameters: the 
mass ratio of PAHs to total dust normalized to the Galactic 
value of 0.046 (/pah): 

1 M PAH o + M PAH+ 



where Mpah° and M PAH + are the masses of the neutral 
and ionized PAHs respectively, and Mp A h = Mpah° + 
M PAH +. 

The stellar continuum makes a large contribution to 
the SED at NIR wavelengths. It needs, as a free parameter, 
the total stellar mass, M i ds t ar - M i datar is used to com- 
pute the stellar emission, L„(A)* 6 . Note, however, that the 
resulting value for Moidstar is not actually a proper measure- 
ment of total stellar mass. The latter would require data at 
shorter wavelengths, different extinction geometries, vary- 
ing the stellar age and eventually the Initial Mass Function 
(IMF). Rather, M i dstar is a parameter that is varied to un- 
derstand the contribution from evolved stars into the MIR. 

The dust must be excited by optical-UV radiation to 
emit. We adopt the shape of the interstellar radiation field 
(ISRF) of the Galactic diffuse ISM by Mathis et al. (1983), 

The adopted ISRF, U, is normalized by the local Solar 
neighbourhood value (taken to be 2.2 x 10~ 5 W m~ 2 ): 

/1(A) dX 



U = 



L 



09/um 



2.2 x 10" 5 



(3) 



where the lower limit on the integral is set at the Lyman 
break and the upper limit approximates the value at which 
excitation is no longer significant. 

The monochromatic luminosity of dust exposed to U 
per unit dust mass, l„(U, A), is: 



lv(U,\) 



dL v (U,X) 
dM d (U) 



(4) 



This monochromatic luminosity can be expressed as a sum 
of the three dust components (PAHs, silicate and graphite): 



l„(U,\)= f P AH 



fpAH^ AH+ {U,\) + {l 



fpAH+) 



x£ AH (U,X) 



+ if(u,\) + C 



(U,X) (5) 



Here l?(\J,\), where N = PAH+, PAH°, sil. and gra., are 
the luminosities of the ionized PAHs, neutral PAHs, silicates, 
and graphites respectively. 

With an admixture of dust components in an interstel- 
lar radiation field, we then need to consider the possible 
range of environments within which these components could 
be found. Following Dale et al. (2001), this can be done via 
the parameter, a, which relates the heated dust mass, M d , 
to the radiation field, U : 



dM d {U) oc U~ a dU 



(6) 



between U m in and Umax- The range of a, which is defined 
by Eqn. 6 is typically 1 to 2.5, depending on whether the 
dust is in a diffuse medium (a w 2.5) or whether the dust 
is in a dense medium (a ~ 1). Since the area measured 
contains contributions from both kinds of media, the a from 



The subscripts v or A indicate that this is a monochromatic 
term, and the A in brackets indicates a wavelength dependence. 
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the model is a characteristic value for the measured region 
and will likely fall somewhere between 1 and 2.5. 

Therefore, with Equations (4) and (6), and the stellar 
component added, the total SED can be expressed as 



L V (X) 



(1 - a)M d r u 

jjl — a _ Tj-l — a I 
Umax u min J 



l„{U,X)U- a dU + L*(A) 

xexp[-r(A)] (7) 



where the optical depth, r = Ay/1.086, (and Ay is the 
visual extinction), and L*(X) comes from Fioc & Rocca- 
Volmerange's (1997) model 7 for a burst of 5 Gyr and a 
mass of stars, M oHstar . The exp [— r(A)] term accounts for 
the extinction in the 'slab' geometry that we have adopted. 
Note that the visual extinction, Ay, is constrained by the 
J, H, and K band ratios. A 'true' value of Ay requires addi- 
tional measurements at visible wavelengths, so this extinc- 
tion should be deemed 'representative', rather than exact, 
and useful for comparing different regions in the galaxy. 
To summarize, the eight parameters the model is fitting 

are: 

• The mass of dust (M4) in units of M©, 

• The ratio of mass of PAHs to mass of dust normalized 
to the Galactic value (0.046) (/pah), 

• The ratio of mass of ionized PAHs to total PAHs 

(fpAH+)i 

• The mass of the old stellar population (M id s tar) in 
units of Mq , 

• Minimum and maximum heating intensities (Umin & 
Umax respectively) normalized to the local radiation envi- 
ronment (2.2 x 10~ 5 W m" 2 ), 

• The visual extinction (Ay), and 

• a defined in Equation 6. 

Finally, we note that the parameters a, Umin and Umax 
do not have a strong physical meaning when taken indepen- 
dently. However, these parameters together quantify the dis- 
tribution of dust mass through different environments. The 
fraction / co id of dust mass exposed to radiation densities 
lower than U < 8.5 (corresponding to a silicate equilibrium 
temperature of 25 K) is: 



/cold — 



U(25 K) 



'l-a)U- 



U, 



1-a 

max 



ui: 



■dU 



(8) 



These quantities will also be calculated from the model. 

This model is integrated within each filter, with the 
proper flux convention, and then compared to the observa- 
tions, by minimizing the \ 2 - In addition, we estimate the 
error on the value of the parameters as follows. For each ob- 
served SED, {F v (X i )} i , with error bars {AF 1/ (X i )} i , we per- 
form several fits of {F v (Xi) + 9i x AF„(Xi)} v where 9i are 
independant normally distributed random variables. There- 
fore, such a fit provides the values of the parameters for one 
particular realization of the error bars. By iterating over 
several 9i sets, we obtain the distribution of each parameter 
value. The dispersion of this distribution provides the error 
bar on the parameter. 



7 The NIR old stellar continuum is roughly independent of the 
age of the populations for a burst > 1 Gyr. The stellar continuum 
is modeled with the population synthesis code PEGASE. 



4.3 SED results for NGC 891 

Figure 12 shows the observed SEDs for the six measured 
regions in NGC 891 as well as the model fits. The yellow line 
represents the stellar contribution corrected for extinction, 
the red line represents the dust and PAH contribution, the 
grey line is the total line of best fit and the blue line is the 
SED without extinction. The data points are circles with 
error bars and the green points are the model results in the 
observed bands. 

In general, the model results agree with the data to 
within estimated errors throughout the NIR to the sub-mm. 
The exceptions are the A 450 /im points, which are often 
too high. One possibility is that the observed discrepancy is 
real and that the dust properties in this part of the spec- 
trum are peculiar. However, it is also possible that the A 450 
/im fluxes may be in error. We have accepted the original 
calibration from previously published results at this wave- 
length (Sect. 2.3) and both Israel et al (1999) and Alton et 
al. (1998), who observed this galaxy with 2 independent sets 
of observations, give consistent flux values. Nevertheless, the 
total error for this point may have been underestimated. 

The FIR flux densities in the models peak at shorter 
wavelengths (hotter) compared to what has been observed 
in many other galaxies (e.g. Dunne & Eales 2001; Bendo et 
al. 2003; Dale et al. 2005; Draine et al. 2007b). For exam- 
ple, a better constrained total SED peaks at around A 100 
/im (Galliano et al. 2008a), whereas many of our SEDs are 
peaking at about A 60 - 70 /im (Figure 12). However, in this 
case, it is the relatively low A 850 /im fluxes that prevent the 
SEDs from having a more powerful cold component. Adding 
additional FIR data in the A 100 - 200 fim range would tend 
to make the SEDs a little colder, but it would also make 
the fit of the SCUBA data even worse. Nevertheless, given 
that each region has data points at the same wavelengths, 
the SED model results for the various regions can still be 
compared to each other in a meaningful way. 

The four in-disk regions all show similarly shaped SEDs 
which peak in the A 70 to 80 /im wavelength regime, includ- 
ing the outer north disk which might have been expected to 
show a 'colder' (longer wavelength) peak given its larger ra- 
dial offset from the center (Table 5). However, as indicated 
in Section 4.1, the position of the outer north disk point is 
96 arcsec from the centre which coincides with a broad peak 
in the Ha emission (see image in Kamphuis et al. 2007). It 
also corresponds to a peak in the surface brightness profile 
along the disk (see Figure 8). Therefore, its SED is likely 
more representative of a SF region than of the far outer re- 
gions of a typical galactic disk. The halo SEDs, on the other 
hand, do peak at longer wavelengths suggesting the domi- 
nance of colder dust at high latitudes. Adopting lower A 70 
/im values should not change this conclusion. 

The halo SEDs differ from the disk SEDs in other sig- 
nificant ways. First of all, they show an overall lower surface 
brightness by about a factor of 10 (see absolute values on 
ordinate axes, the apertures being indentical), as might be 
expected for a faint halo region. They also reveal that virtu- 
ally no extinction (by our representative Av, see Section 4.2) 
is required to match the NIR emission which is dominated 
by starlight. Again, this is not surprising, given the contrast 
between observing stellar emission through a dense dusty, 
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edge-on disk, or observing starlight in a faint halo region 
with mainly diffuse tenuous dust. 

The other contrast between the halo and disk SEDs 
are in the relative brightness of the MIR in comparison to 
other parts of the spectrum. The halo SEDs show a lower 
brightness in the MIR, where PAHs are the dominant con- 
tributor, in comparison to disk SEDs. For example, the three 
in-disk regions (not including centre) all show similar ratios 
of PAHs to old stars, as measured by the A 7.7 fim peak to 
the dust-corrected A 1 /im peak ratio, ranging from ~ 0.6 to 
0.9. However, this ratio in the halo is somewhat depressed 
in comparison to the disk, with values in the range w 0.2 to 
0.4. As we showed in Figure 10, the stellar halo, as measured 
by the K-band flux, is also broader than the A 7.7 /im halo 
emission. Therefore, the lower PAH/stellar ratio in the halo 
is consistent with the vertical distribution seen earlier. 

PAHs in the halo are also depressed in comparison to 
cold dust. The three in-disk regions (again, not including 
the central region) show similar PAH/cold dust ratios, as 
measured by the A 7.7 /im peak to the FIR peak ratio. We 
cannot compare this ratio directly for the halo SEDs since 
the latter do not have well defined peaks, with only upper 
limits at A 70 fim. Using the A 7.7 jj,m/X 850 /im ratio instead 
to probe the PAH/cold dust emission, we find that the disk 
SEDs have significantly higher ratios, of order 2000, whereas 
the halo ratios are of order 300. 

Thus, there is a deficit of PAH emission in the halo 
with respect to the stars and cold dust, as our global halo 
profiles from Sect. 3.3 suggested (Figures 10 and 11). In 
fact, the PAH strength is also lower in the halo than in the 
centre of the galaxy, though this is not quite as obvious from 
Figure 12. Rather, it is a result of the model fitting, the full 
results of which we provide in Table 6. 

From Table 6 we find that there is no significant dif- 
ference between M<j, / P/U} +, a, Av and Moidstar for the 
disk regions, although M id s tar is predictably higher for the 
central region. Except for a, these all of the parameters are 
also predictably lower for the halo regions. However, /pah, 
U m in, and U m ax, and f co id show significant variation, in gen- 
eral. 

We will discuss all of these results more fully in Sects. 
5.1 and 5.2. 

From the SED model, the masses of the PAHs, Mpah, 
for the six regions are implicitly included in Table 6 from 
Mpah = 0.046 fpAH M d . As with the other modeled 
masses given in Table 6, PAH masses will be specific to the 
adopted size and location of the region. Nevertheless, it is 
significant that these quantities can be calculated from the 
model. Few mass estimates for PAHs in external galaxies 
exist at this time. Indeed, we are the first to determine PAH 
mass estimates specifically for a halo region. 
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Figure 12. SEDs for the various regions of NGC 891 given in Table 5, fit with the model described in Sect. 4. The yellow line represents 
the stellar contribution corrected for extinction, the red line represents the dust and PAH contribution, the grey line is the total line of 
best fit and the blue line is the SED without extinction. The data points are circles with error bars and the green points are the model 
results in the observed bands. 



Table 6. Modeling Parameter Results for Different Regions of NGC 891. 



parameter 01 '' 


centre (1) 


south disk (2) 


inner north disk (3) 


outer north disk (4) 


M d (Mq) 
fpAH 
fpAH+ 
a 

U Min 
Um ax 

A v 

M l ds tar (Mq) 

fcold (percentage) 


(7.3 ± 1.9) x 10 5 (26%) 
0.9 ± 0.2 (25%) 

nfi-Ln 1 ('\ ^9?,^ 

2.2 ± 0.5 (24%) 

5 ± 2 (36%) 

654^f (285%) 

5.7 ± 0.2 (3%) 

(1.1 ±0.0) x 10 10 (3%) 

51 ± 25 (48%) 


(4.9 ± 1.4) X 10 5 (28%) 
1.5 ±0.4 (28%) 
ne-Ln 1 ('\ 

2.4 ± 0.2 (8%) 

3 ± 1 (32%) 

36860 ± 30208 (82%) 

4.0 ±0.1 (3%) 

(3.8 ±0.1) x 10 9 (3%) 

80 ± 11 (14%) 


(5.3 ± 1.3) x 10 5 (25%) 

1.8 ± 0.4 (25%) 

U.D I U.l yZL /§) 

2.4 ±0.1 (4%) 

2 ± 1 (29%) 

21481 ± 23816 (111%) 

3.9 ±0.1 (3%) 

(3.8 ±0.1) x 10 9 (3%) 
86 ± 6 (7%) 


(7.8 ± 3.9) x 10 6 (51%) 

1.8 ± 0.4 (21%) 
u.o zn u.i {zu/o) 

1.9 ± 0.2 (13%) 
± 1 (171%) 
601 ± 521 (87%) 
4.3 ± 0.2 (4%) 

(2.8 ±0.1) x 10 9 (4%) 
99 ± 1 (1%) 


parameter 


inner halo (5) 


outer halo (6) 






M d (Mq) 
fpAH 
fpAH+ 
a 

UMin 
UMax 

A v 

Moldstar (Mq) 

fcold (percentage) 


(2.5 ± 0.9) X 10 5 (35%) 

1.5 ±1.7 (117%) 

0.4 ± 0.3 (76%) 

2.5 ±0.1 (2%) 

2 ± 2 (110%) 

18914 ± 32890 (174%) 

0.0 ± 0.0 (422%) 

(9.6 ± 0.2) x 10 8 (2%) 

93±I 6 (17%) 


(1.5 ±0.7) x 10 5 (45%) 

1.2 ± 1.5 (120%) 
0.9 ± 0.2 (18%) 

2.3 ± 0.2 (8%) 
2 ± 3 (106%) 
54711 ± 47703 (87%) 
0.0 ±0.1 (148%) 

(4.2 ±0.1) x 10 8 (2%) 
84±£ (26%) 







a See text Sect. 4.2 for the definition of the parameters. 

^ The parameters a, UMin an d l?Max are correlated. That is the reason why they span a wide range of values and are not normally distributed. 
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5 DISCUSSION 

5.1 The Dusty Disk of NGC 891 

The changing morphology of the disk emission of NGC 891 
as a function of wavelength, as indicated in Sect. 3.1, nicely 
illustrates the changing admixture of IR-emitting compo- 
nents — from dominant starlight at NIR wavelengths to 
PAHs and dust as the wavelength increases. The PAH/dust 
distribution is complex, as shown in the maps (Figures 2 - 
4, Al - A4) as well as the high resolution, high S/N A 8 /j,m 
PAH major axis profile (Figure 8) . The latter figure reveals 
new details of the PAH distribution along the major axis. 

It is useful to compare the major axis profile of Figure 
8 to the molecular gas distribution. Scoville et al. (1993) 
showed that the molecular gas in NGC 891 has a peak at 
the center surrounded by a molecular ring located 3-7 kpc 
(65 - 150 arcsec) in radius, a structure also detected by Sofuc 

6 Nakai (1993). Our PAH profile, which reveals much more 
detail, follows this general molecular gas distribution quite 
well, showing a central peak and two peaks on either side at 
similar radii. In addition, we detect an A 8 /im peak on the 
south-west side of the disk at a radius of « 140 arcsec (6.5 
kpc) which was seen in CO by Sofue & Nakai (1993). Thus, 
the dominant correlation of the PAH emission in NGC 891 
is with the molecular gas, a result that is in agreement with 
Le Coupanec et al. (1999) and also seen in other galaxies 
(e.g. Irwin & Madden 2006; Regan et al. 2006). The much 
lower PAH emission at larger radii is likely due to lower 
quantities in the lower density extra-planar gas as well as 
less excitation farther from SF regions. 

Although our data were not of sufficient quality to ob- 
tain SEDs at every point in the disk, we have obtained good 
results for four locations as listed in Table 6. 

The ratio, fpAH, is given in terms of the MW value 
of 0.046. Thus, we can immediately compare the regions in 
NGC 891 to the MW in this respect. For the four in-plane 
regions of NGC 891 this ratio ranges from about 0.9 — > 2, 
indicating that the PAH abundance is roughly Galactic. Al- 
though the 'centre' of this edge-on galaxy also includes com- 
ponents along the line of sight, SED results for the centre 
should nevertheless be weighted towards the true centre of 
the galaxy and the fact that the central SED differs from the 
other SEDs supports this view. For example, we find that 
fpAH at the centre of the galaxy is the lowest of the four in- 
plane regions. This agrees with the results of Le Coupanec et 
al (1999) who found a reduction in the A 11.3 fim PAH fea- 
ture at the galactic centre. It is also consistent with Smith et 
al. (2007) and Bendo et al. (2008) who both find a decrease 
in PAH emission in the centres of some galaxies. The fraction 
of dust exposed to low ISRFs (/ co id) is — 50% in the centre, 
where a significant fraction of the dust emission comes from 
hotter dust, likely embedded in star forming regions / co id is 
between ~ 80% and 100% in the disk, where the emission 
seems to come primarily from a diffuse component. 

Although the strength of the ISRF at lower wave- 
lengths, Umin is approximately the same for all in-disk re- 
gions, surprisingly, the heating intensity, U ma x, is greatest 
in the south disk, suggesting a greater range of heating en- 
vironments than in the other in-disk locations. As indicated 
in Sect. 4.3, there is an enhancement in H a emission at this 
location and the higher Umax likely reflects the fact that the 
south disk samples a region of star formation. These results 



also suggest that the enhanced H a emission is likely real 
and not only the result of dust obscuration elsewhere (see 
discussion in Kamphuis et al. 2007). Note however, that the 
Umin and Umax values are constrained only by the four data 
points in the FIR, and therefore may reflect some of the 
noise in those values. 



5.2 Extraplanar Dust and PAHs in NGC 891 

Previous investigations have shown that NGC 891 has ex- 
traplanar emission in neutral gas, DIG, molecular gas, and 
dust (eg. Oosterloo et al. 2007, Fraternali et al. 2005, Rossa 
et al. 2004, Garcia-Burillo et al. 1992, Howk & Savage 2000, 
see Sect. 1 and Table 7). Our results now show that PAHs 
can be added to this list, with a measured z height extent of 
at least 2.3 kpc (Sect. 3.3). The z width of the cold dust dis- 
tribution at A 450 /im is greater than those of both the warm 
dust (A 24 lira) as well as PAHs (A 7.7 and 8 /im), as shown 
in Figure 10 and noting error bars (see also Figure 11). 

Warm dust is only seen close to heating sources and 
therefore it is not surprising that this component has a nar- 
rower z width than cold dust. As for PAHs, we expect that 
PAHs may very well exist at even higher z latitudes than ob- 
served because we know that PAHs are present in H I (see 
the SED models of Dwek et al. 1997, Draine & Li 2007a, and 
Zubko et al. 2004 which show prominent PAH features asso- 
ciated with the diffuse ISM) and the H I halo of NGC 891 is 
extensive (Table 7) . The fact that the PAH vertical width is 
closer to that of the warm dust (Figure 10) therefore suggests 
that the dominant PAH excitation sources are also embed- 
ded in the disk. Recent Galaxy Evolution Explorer obser- 
vations have now revealed an FUV halo in NGC 891 (Gil 
de Paz et al. 2007) which potentially provides a halo source 
of FUV photons. However, the FUV emission in galaxy ha- 
los is largely due in-disk stellar continuum photons that are 
reflected from halo dust 8 (Hoopes et al. 2005). In an op- 
tically thin halo environment, FUV photons from the disk 
should dominate over scattered photons and, in either case, 
the dominant excitation source ultimately resides from stars 
the disk. 

The extents of the extraplanar emission are summarized 
in Table 7. These components are listed in order of largest 
distance from mid-plane to smallest distance. All of the val- 
ues given are for the farthest measured extent of the halo 
(as opposed to a characteristic scale height) and thus, they 
depend on whatever S/N is available for the data. Therefore, 
Table 7 is meant to provide a measure of the full extent of 
the extraplanar components in NGC 891 to the limits of cur- 
rent data. Nevertheless, Figures 10 and 11 display how the z 
width of a subset of these quantities compare to each other 
and the widths of these profiles agrees with the ordering of 
Table 7. 

An important result of this work is our first attempt at 
a SED model in the extraplanar (halo) region of an exter- 
nal galaxy (Table 6). The extra-planar dust mass results are 
consistent with those determined from extra-planar absorp- 
tion (Howk & Savage 1997), after reasonable corrections for 
different apertures are taken into account. The /pah ratios 
in the halo are similar to those in the disk. This result was 

8 There could also be a minor contribution from shock excitation. 
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Table 7. Composition and Extent of Extraplanar Components 
in NGC 891. 



ISM component distance from reference 
mid-plane (kpc) 

neutral gas 22 Oosterloo et al. (2007) 

15 Fratcrnali ct al. (2005) 

5 Swaters et al. (1997) 

old stars 3 this paper (Figure 10) 

DIG 2.2 Rossa ct al. (2004) 

2 Howk & Savage (2000) 

warm dust 2.5 this paper (Figure 10) 

cold dust 2 Alton ct al. (1998) 

PAHs 2.3 this paper (Figure 10) 

molecular gas 1 - 1.4 Garcia-Burillo et al. (1992) 



gas, as measured by A 450 /jm emission, is broader in z than 
both warm halo gas A 24 (im and PAH emission, suggesting 
that the dust heating sources and PAH excitation sources 
are within the disk. 

We have used our multi- frequency data set to construct 
complete IR spectra for four different locations in the disk 
of NGC 891 and two in its halo. This has allowed us to 
model SEDs in different environments both within the disk 
of NGC 891 and, for the first time, in the halo. The SED 
models have allowed us to determine the mass of dust, frac- 
tion of PAHs, and fraction of ionized PAHs for these differ- 
ent environments (Table 6). We find that the PAH fraction, 
fpAH, is similar to Galactic values (within a factor of two) 
and this value is lowest in the galaxy's center, consistent 
with previous studies of PAHs in galaxies. The fraction of 
dust exposed to a colder (i.e. equivalent to an equilibrium 
dust temperature less than 25 K) radiation field, f co id, is 
also lowest in the center, suggesting the presence of a hotter 
radiation field at the center of the galaxy. 

In the extraplanar dust, /pah is similar to fpAH in the 
disk. Therefore we cannot conclude a lower PAH abundance, 
especially since the SED model used does not take into ac- 
count a changing radiation field that is more typical of a 
halo. The model used could take into account the changing 
radiation field, but we would have needed additional con- 
straints. The fraction of PAHs exposed to a cold radiation 
field is very high in the halo and is almost unity. Together 
with the vertical halo extents indicated above, these results 
suggest that the halo PAHs are excited from escaping disk 
photons. 
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unexpected given the observations of the vertical profiles in 
Figure 10, which show that the cool dust (at A 450 /im, for 
example) extends further into the halo vertically than the 
PAHs (at A 7.7 fim, for example). However, the error bars 
are very large for these halo results, and the SED model 
does not take into account the variation of the stellar popu- 
lations from the disk to the halo. A decrease of the hardness 
of the ISRF outward, would decrease the PAH strength, for 
a constant PAH-to-dust mass ratio, which is what we see 
here. The quantity, f co id, which is close to 1 in the halo, is 
consistent with an interpretation in which few hot photons 
are available for excitation. 

While NGC 891 is known for large amounts of neutral 
gas, DIG, and dust in its halo (Sect. 1), its PAH emission 
is not as extended in z as has been seen in several other 
galaxies which are not known for abundant halo gas and 
dust (see Irwin & Madden 2006, Irwin et al. 2007). Since 
the excitation source for the PAHs appears to be within the 
disk (see above), we speculate that FUV photons that could 
potentially excite high latitude PAH emission in NGC 891 
may instead be absorbed by other gaseous and dust compo- 
nents. That is, other components in the halo of NGC 891 
have some optical depth to FUV photons that originate in 
the disk. 



6 CONCLUSIONS 

We have examined 20 spatially resolved IR data sets for 
NGC 891, 14 of which are newly reduced and/or newly pre- 
sented data. The emission spans wavelengths from A 1.2 to 
850 /jm. Although there are a variety of resolutions and S/N 
for the different data sets, which include both ground-based 
and space-based observations, we have compared the spa- 
tial distributions, both along and perpendicular to the ma- 
jor axis for a subset of the data. PAH emission is found to 
correlate well in the disk with molecular gas, likely because 
of higher densities of gas and dust in these regions as well 
as the abundant supply of UV photons from SF regions. 

Of particular interest is the extraplanar, or 'halo' gas 
and dust in NGC 891. The H I, DIG, and dusty halos have 
been established by others. Our work has compared the z 
extents of various components and also shown that PAHs 
also exist in the halo of this galaxy to heights of 2.3 ± 0.2 
kpc at a 3cr level. NGC 891 is only the 4th galaxy in which 
PAHs have been detected in the halo. We find that cool halo 
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Figure Al. Surface brightness contour maps for remaining ISO bands that isolate PAH emission, where the others were shown in Figure 
2. Symbols as in Figure 2. Contours arc 0.3 (3a), 0.8, 2, 4, 6, and 8 MJy/sr for A 6.0, ^m and 0.7 (3<r), 1.5, 4, 6, 8, 10 MJy/sr for A 11.3 
fim. 
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Figure A2. Surface brightness contour maps for the remaining ISO and IRAC wave bands that isolates stellar emission, where the other 
was shown in Figure 3. Symbols as in Figure 2. Contours are 0.08 (3<r), 0.15, 0.4, 1.2, 2.5, and 5 MJy/sr for A 3.6 and 0.3 (3a-), 0.7, 
1.5, 3, 4, 6 MJy/sr for A 4.5 fim. 
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Figure A3. Surface brightness contour maps for the remaining ISO wave bands that isolate the MIR continuum, where the other was 
shown in Figure 4. Symbols as in Figure 2. Contours arc 0.5 (3<r), 1.5, 4, 6, 8 MJy/sr for A 6.8 jum, 0.3 (3<r), 0.8, 2, 4, 6, 8 MJy/sr for 
A 9.6 Atm, and 0.5 (3ct), 1.5, 4, 6, 8 MJy/sr for A 14.9 ^tm. 
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Figure A4. Surface brightness contour maps for the ISO wide bands. Symbols as in Figure 2. Contours are 0.3 (3cr), 0.8, 2, 5, 8 MJy/sr 
for A 6.7 fim and 0.5 (3a), 1, 2, 4, 6, 8 MJy/sr for A 14.3 ^tm. 
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Figure A5. The SCUBA A 450 fim and A 850 fim maps showing cool dust in NGC 891 (see Alton et al. 1998). Symbols are as in Figure 
2. For the A 450 /an map, the contours are 2.5, 4, 6, 8, 9.5 Jy beam -1 . For the A 850 £tm map, the contours are at 31 (3<r), 50, 80, 125, 
190, and 250 mjy beam -1 . 



